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Abstract — For the nonlinear power flow problem specified with 
standard PQ, PV, and slack bus equality constraints, we present 
a sufficient condition under which the specified set of nonlinear 
algebraic equations has no solution. This sufficient condition is 
constructed in a framework of an associated feasible, convex op- 
timization problem. The objective employed in this optimization 
problem yields a measure of distance (in a parameter set) to the 
power flow solution boundary. In practical terms, this distance is 
closely related to quantities that previous authors have proposed 
as voltage stability margins. A typical margin is expressed in 
terms of the parameters of system loading (injected powers); 
here we additionally introduce a new margin in terms of the 
parameters of regulated bus voltages. 

Index Terms — Power flow. Power flow solution existence. Max- 
imum loadabiflty, Solution space boundary 



I. Introduction 

POWER flow studies are the cornerstone of power system 
analysis and design. They are used in planning, operation, 
economic scheduling, transient stability, and contingency stud- 
ies 12 • The power flow equations used in these studies model 
the relationship between voltages and active and reactive 
power injections in a power system. The nonlinear power 
flow equations may not have any solutions (the power flow 
equations are said to be insolvable). That is, it is possible 
to choose a set of power injections for which no valid 
corresponding voltage profile exists. Practical cases that may 
fail to have a solution include long-range planning studies 
in which the studied system may not be able to support 
projected loads, and contingency studies for which the loss of 
one or more components may yield a network configuration 
that is similarly inoperable for the specified injections. This 
paper presents a practically computable sufficient condition, 
that when satisfied, rigorously classifies a specified case as 
unsolvable. This method also provides controlled voltage and 
power injection margins that characterize a distance to the 
power flow solvability boundary. 

In engineering practice, large scale nonlinear power flow 
equations are typically solved using iterative numerical tech- 
niques, most commonly Newton-Raphson or its variants [2l, 
that rely on an initial guess of the solution voltage magnitudes 
and angles. An important limitation of these techniques is 
that they are only locally convergent. That is, they generally 
do not converge to a solution from an arbitrary initial guess 
|[T), and may show very high sensitivity and highly complex 
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behavior with respect to initial conditions for certain study 
cases. It is well recognized that a power flow problem may in 
general have a very large number of solutions; for example, 
the work of [3J establishes cases for which the number of 
solutions grows faster than polynomial with respect to network 
size. For cases having multiple solutions, each solution has 
a set of initial conditions that converges to that solution 
in Newton-Raphson iteration. Characterization of Newton- 
Raphson regions of attraction was the subject of |4|, which 
demonstrated cases for which the boundaries of these attractive 
sets were factual in nature. So despite the fact that very large 
scale problems (lO's or lOO's of thousands of unknowns) are 
routinely solved in power engineering practice, it is important 
to recognize that as parameters move outside of routine 
operating ranges, the behavior of these equations can be highly 
complex, and failure of convergence for a Newton-Raphson- 
based commercial software package is far from a reliable 
indication that no solution exists. 

The properties of the Newton-Raphson iteration guarantee 
(under suitable differentiability assumptions) that the iteration 
must converge to the solution for an initial condition selected 
in a sufficiently small neighborhood about that solution [5|. 
However, when a selected initial condition (or some limited 
set of multiple initial conditions) fails to yield convergence, 
the user of a Newton-Raphson-based software package is left 
with an indeterminate outcome: does the specified problem 
have no solution, or has the initial condition(s) simply failed 
to fall within the attractive set of a solution that does exist? 

Conditions to guarantee existence of solutions to the power 
flow equations has been an active topic of study. For example, 
|6| describes sufficient conditions for power flow solution 
existence. However, as sufficient conditions, these are often 
conservative: a solution may exist for a much larger range 
of operating points than satisfy the sufficient conditions. 
Other work on sufficient conditions for power flow solvability 
includes |7|, which focuses on the decoupled (active power- 
voltage angle, reactive power-voltage magnitude) power flow 
model. Reference [8] describes a modified Newton-Raphson 
iteration tailored to the type of ill-conditioning that can 
appear in power systems problems. While convergence to a 
solution may be judged a constructive sufficient condition to 
demonstrate solvability, such approaches do not escape the 
fundamental limitations of a locally convergent iteration. In 
more recent work, [9] provides two necessary conditions for 
saddle-node bifurcation based on lines reaching their static 
transfer stability limits; however, this work does not yet 
provide a test for power flow solvability or define a distance 
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to the power flow solvability boundary. 

A measure of the distance to the solvability boundary (the 
set of operating points where a solution exists, but small 
perturbations may result in the insolvability of the power flow 
equations fTU\) is desirable to ensure that power systems are 
operated with security margins. If a solution does not exist for 
a specified set of power injections, a measure of the distance 
to the solvability boundary indicates how close the power flow 
equations are to having a solution. Existing work in this area 
uses a Newton-Raphson optimal multiplier approach [11] to 
find the voltage profile that yields the closest power injections 
to those specified [12 1, |13|. Techniques for finding load 
margins comprised of the minimum distance to the power flow 
solvability boundary include [14] and [Tsl. An algorithm that 
combines continuation and non-linear optimization techniques 
to either solve the power flow equations, when possible, or 
calculate a measure of power flow insolvability is presented 
in [16|. Reference |17| describes an optimization problem 
that applies interior point methods to minimize the load 
shedding necessary to obtain solvable power flow equations. 
The minimum amount of load shedding is used as a measure 
of power flow insolvability. Investigating the worst-case load 
shedding necessary for power flow solvability is also discussed 
in |18| and [19]. Reference [20| summarizes and compares 
some of these power flow insolvability measures. 

In this paper, we present a sufficient condition under which 
a specified power flow problem is guaranteed to have no 
solution. The computation yields as a by-product controlled 
voltage and power injection margins to the power flow solv- 
ability boundary. 

The sufficient condition for power flow insolvability is 
based on an optimization problem that includes a relaxation 
of certain equality constraints in the power flow equations. 
Specifically, in this optimization problem, the voltages at slack 
and PV buses are not fixed, but instead have a one-dimensional 
degree of freedom (i.e., they are allowed to change in constant 
proportion). In Section |III] we prove that the extra degree of 
freedom guarantees the modified power flow equations to have 
at least one solution. In an idealized lossless case, one may 
interpret this as follows: a sufficiently high voltage profile 
allows the system to meet any specified power injections. 
By continuity from the lossless case, we argue that this will 
continue to hold for modest losses, as is typical of models 
for bulk transmission. With the relaxed problem feasible for 
some (sufficiently high) voltage profile, we establish a non- 
empty feasible set for the optimization. 

With a non-empty feasible set established, the optimization 
problem then seeks to minimize the slack bus voltage mag- 
nitude (using the one-degree-of-freedom in voltage profiles), 
subject to the active and reactive power injection constraints 
of the power flow equations. Importantly, we will show 
that a further relaxed version of this optimization problem 
is a convex semidefinite programming problem, and hence 
has a practically computable global minimum. If the global 
minimum slack bus voltage obtained from this optimization 
problem is greater than the originally specified slack bus 
voltage, there can be no solution to the originally specified 
power flow equations. However, due to the nature of the 



relaxation, one may not draw a firm conclusion from the 
converse: if the minimum slack bus voltage is less than or 
equal to the specified slack bus voltage, the power flow 
equations may or may not be solvable. 

The ratio of the specified slack bus voltage to the minimum 
slack bus voltage gives a "controlled voltage margin" to the 
power flow solvability boundary. In a provably unsolvable 
case, this margin is the multiplicative factor by which the 
controlled voltages must be increased to allow the possibility 
of existence for a power flow solution. 

It is widely recognized that the power flow equations are 
quadratic in the complex voltage vector when these voltages 
are expressed in rectangular form. Exploiting this fact, an 
analogous power injection margin can also be calculated; 
here the new, one degree of freedom introduced represents 
a constant power factor change in injections at each bus in 
proportion to the specified injections. When the power flow 
equations do not have a solution, the power injection margin 
provides the factor by which the power injections must be 
decreased to admit the possibility for power flow solution 
existence. 

These margins are non-conservative upper bounds. Thus, 
for an insolvable set of specified values, a change in voltage 
by at least the amount indicated by the voltage margin (or a 
change in power injections by at least the amount indicated 
by the power injection margin) is required for the power flow 
equations to be potentially solvable. 

The dual of the optimization problem used in the sufficient 
condition can be written as a semidefinite program (SDP). The 
optimal power flow problem (i.e., finding the optimal operating 
point for a power system subject to physical and engineering 
constraints) was recently formulated as an SDP (2V\, f22]. 
In prior work, the authors created an SDP formulation of 
the power flow equations in an attempt to calculate multiple 
solutions to these equations f23\. In contrast to the non-convex 
primal optimization problem [24), the feasible region of the 
dual problem formulated as an SDP is convex. The optimal 
objective value obtained from the dual SDP formulation is 
a lower bound on the objective function value used in the 
sufficient condition here. Thus, if the sufficient condition holds 
based on the lower bound from the dual SDP formulation, one 
can be assured the originally formulated power flow equations 
admit no solution. 

The organization of this paper is as follows. In Section 
Hn we give an overview of the power flow equations. In 
Section |III1 we provide the existence proof that shows the 
feasibility of the optimization problem used by the proposed 
condition. In Section IIVI we describe the sufficient condition 
for power flow insolvability and define voltage and power 
injection margins. Numeric examples are provided in Section 
rv] Section |VT] examines solutions with non-zero duality gap. 
We then conclude with a discussion of future work. 

II. Power Flow Equations Overview 

The power flow equations describe the sinusoidal steady 
state equilibrium of a power network, and hence are formu- 
lated in terms of complex "phasor" representation of circuit 
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quantities (see, for example, Ch. 9 of (25]). The underlying 
voltage-to-current relationships of the network are linear, but 
the nature of equipment in a power system is such that 
injected/demanded complex power at a bus (node) is typically 
specified, rather than current. The relation of interest is be- 
tween the active and reactive power injected at each bus and 
the complex voltages at each bus, and hence the associated 
equations are nonlinear. Using the standard polar represen- 
tation for complex voltages, and rectangular "active/reactive" 
representation of complex power, the power balance equations 
at bus i are given by 



P^^V,}_^ Vk (G,fc cos (S, - 4 

fc=i 



fc=i 



- B,k sin (Si ~ Sk)) (1) 



Bik cos {S^ - Sk)) (2) 



where Pi and Qi are the active and reactive power injections, 
respectively, at bus i, Y — G + jB is the network admittance 
matrix, and n is the number of buses in the system. 

To represent typical behavior of equipment in the power 
system, each bus is classified as PQ, PV, or slack, according to 
the constraints imposed at that bus. PQ buses, which typically 
correspond to loads, treat Pi and Qi as specified quantities, 
and enforce the active power ([TJ and reactive power (|2} 
equations at that bus. PV buses, which typically correspond to 
generators, specify a known voltage magnitude Vi and active 
power injection Pi, and enforce only the active power equation 
(HJ. The associated reactive power Qi may be computed as an 
"output quantity," via (|2]). Finally, a single slack bus is selected, 
with its specified Vi and 6i (typically chosen to be 0°). The 
active power Pi and reactive power Qi at the slack bus are 
determined from ([T) and (|2); network-wide conservation of 
complex power is thereby satisfied. 

III. Solution Existence Proof 

The sufficient condition for power flow insolvability re- 
quires the evaluation of an optimization problem in which 
the feasible set is defined by a modified form of the power 
flow equations. The modification introduces one new degree 
of freedom, allowing voltage magnitudes at the slack and PV 
buses to vary; this variation is restricted to a one-degree- 
of-freedom "ray," with all voltage magnitudes changing in 
constant proportion to their base-case values. We prove that 
the feasible space is non-empty for any lossless power system 
(i.e., a network model in which line conductances are all zero). 
Using standard results of basic circuit theory and continuity, 
we argue that the problem retains a non-empty feasible set 
when perturbed with small line conductances, as are typical 
in bulk transmission. 

The proof of solution existence may be outlined as follows. 
We first establish that a solution must exist for any lossless 
system with zero power injections. We then use the implicit 
function theorem to establish that solutions continue to exist 
for injections within small ball around zero. Hence, within this 
ball must exist a ray that aligns with the originally specified 
vector of non-zero power injections. We exploit the quadratic 



nature of the power flow equations to "scale up" voltage 
magnitudes along our one-degree-of-freedom, observing that 
the power injections must likewise move along the previously 
identified ray. It follows trivially that there will exist a scaling 
of voltages such that the specified power injections are real- 
ized, thereby yielding a solution to our modified form of the 
power flow equations. 

A. Existence of a Zero Power Injection Solution 

Consider a generic lossless power system with all active 
and reactive power injections at PQ buses set to zero and all 
active power injections at PV buses set to zero. As our goal 
is accomplished if we can establish existence of one solution, 
we restrict attention to candidate solutions in which all buses 
have the same voltage angle of zero. 

First, since zero power injection at a PQ bus implies 
zero nodal current injection, such buses have only branch 
admittances incident (i.e., from a circuit perspective, these 
are nodes with no independent source connected to them). 
They can be eliminated from the network, and the network 
admittance matrix algebraically reduced via standard results of 
linear circuit theory. We generically assume that the reduced 
network does not result in any zero impedance linesQ 

Next, the substitution theorem [25] guarantees that at any 
PV bus that has an associated non-zero reactive power injec- 
tion, there must exist a shunt admittance of appropriate value 
such that, when substituted in place of the reactive injection, 
an identical solution for bus voltages is preserved. The fact 
that the injections being replaced are purely reactive ensures 
that the associated admittances will be purely imaginary; i.e., 
susceptances only. 

With PQ buses eliminated and reactive injections at PV 
buses replaced by equivalent susceptances, the resulting net- 
work has the property that active and reactive power injections 
at all non-slack buses are identically zero. This allows us to 
write the remaining network constraints of interest as linear 
voltage/current relationships, as follows 



-^slack 






jb2 









jBs +idiag(Ad) 



Vpv 



(3) 



where Ad is is a vector of shunt element susceptances, 
diag (Ad) denotes the diagonal matrix with elements of Ad 



on the diagonal, B 





b2 ' 




B3 



is the bus susceptance 



matrix, and superscript T indicates the transpose operator 
Vsiack and Isiack are the voltage and current injection at the 
slack bus, respectively, and Vpv is the vector of PV bus 
voltages. Note that the lossless assumption implies that the 
network admittance matrix is purely imaginary. 
Solving (O for Ad yields 



Ad = (diag (ypv))"' (-fo2l4iack - Bal/pv) 



(4) 



Because the voltage profile solution we seek is restricted 
to have the same voltage angle at all buses and a non-zero 

'Such a zero impedance line outcome can be eliminated by an arbitrarily 
small perturbation to the underlying line parameter data. 
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voltage magnitude at the slack bus, it follows that the voltage 
at every bus must be non-zero and diag (Vpv) is invertible. 
Hence, for a lossless system under the assumptions specified, 
Q yields a unique solution for the shunt susceptance values 
whose existence follows from the Substitution Theorem. 

provides a zero power injection 



Thus, the vector 



solution to the reduced network that resulted from elimination 
of PQ buses; voltages at PQ buses can be trivially recon- 
structed. We conclude that any lossless system is guaranteed 
to have a zero power injection solution. 

To illustrate that this need not be the case for systems with 
large conductive elements in their bus admittance matrix (i.e., 
high transmission losses), consider the two-bus system with a 
slack bus and a PV bus shown in Fig. [T] 



VsiackZ.Odeg 



jb 



Vpv Z3 deg 



Fig. 1. Two-Bus System 

The transmission line admittance is g + jb; note that in 
this admittance representation, the conductive term g and 
the susceptance term jh appear as parallel branch elements 
between the two buses. The voltage at the slack bus is denoted 
by Vsiack, and the voltage at the PV bus is represented by Vpv 
with angle 9. The power injection at the PV bus is 



PV 



ckigcosi9)+bsm{0)) (5) 



The two-bus system has a zero power injection solution for 
a given set of parameters g, b, Vpv, and Vsiack if a value of 
6*0 exists such that Ppv (do) = 0. The existence of such a 
value of Of) depends on the ratio of Vpv to Vsiack and the 
ratio of b to g. A zero power injection solution to this system 
exists when line resistances are small relative to line reactances 
and voltage magnitude differences are small. Specifically, the 
existence of a zero power injection solution for the system in 
Fig. [U requires 



Vpv 



lack 



< 1 



(6) 



Since voltage magnitudes differences and line resistance to 
reactance ratios are typically small in realistic power systems, 
we expect that typical systems must have zero power injection 
solutions. Consistent with this observation, all the IEEE power 
flow test cases of |26| have zero power injection solutions. 
However, (|6]l confirms that the two-bus example will fail to 
have a zero injection solution when the conductance values 
relative to susceptance is sufficiently large; a zero injection 
solution fails to exist in the case of transmission parameters 
associated with very high losses. 



B. Implicit Function Theorem 

We next apply the implicit function theorem fTf\ at the zero 
power injection solution. Application of the implicit function 
theorem requires a non-singular Jacobian at the zero power 
injection solution. We therefore first investigate the Jacobian 
evaluated at this solution. 

The Jacobian obtained using polar voltage coordinates at 
a zero power injection solution of a lossless system can be 
written as 



Jii 
J22 



(7) 



where 



"11 



J22 — T^TT 



dP 

'85 
dQ 

dV 



-diag (V) (B diag {V) - diag (BV)) (8) 
-diag (V^)B- diag (By) (9) 



Since active and reactive power injections at the slack bus 
are unconstrained, the rows and columns corresponding to the 
slack bus are removed from both Jn and J22. Similarly, since 
the reactive power injections at PV buses are unconstrained, 
the rows and columns corresponding to PV buses are removed 
from J22. Note that both ^ and ^ equal zero for the voltage 
profile with the same voltage angle at all buses corresponding 
to a zero power injection solution of a lossless system. 

The implicit function theorem can be applied at a zero 
power injection solution so long as the Jacobian at this solution 
is non-singular In the lossless case, this requires that J in 
^ is non-singular. We next show that the Jacobian for a 
lossless power system is non-singular at a zero power injection 
solution, provided that all lines are inductive and that the 
network is connected (i.e., no islands). 

J in d?) is non-singular if both Jn and J22 are non-singular 
The matrix diag (BV^) in J22 is equivalent to diag (/), where / 
is the vector of current injections. Since all rows and columns 
in J22 correspond to PQ buses with zero current injections, 
this term is zero. The matrix diag (V) is non-singular since 
all voltages are non-zero for the voltage profile with the same 
voltage angle at all buses. With the slack bus row and column 
removed, B is non-singular for a connected power system with 
inductive lines. Thus, J22 is non-singular 

Since diag {V) is non-singular, Jn is non-singular if 



(Bdiag(V) - diag(BV)) = 
612^2 + ... + Bi„K ••• 



B„iVi 



B„(„_i) V(„_i 



(10) 



is non-singular Note that the diagonal elements in (fTOl i are 
the negative of the sum of the off-diagonal elements in the 
corresponding row. Under the assumption of inductive lines 
such that all off-diagonal elements are negative, this matrix 
has weak diagonal dominance. With the slack bus row and 
column removed, the remaining matrix has at least one row 
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where the diagonal element is strictly greater than the sum 
of the off-diagonal elements (i.e., strict diagonal dominance 
exists for this row). Since the power system associated with the 
B matrix is connected (i.e., no islands), the digraph associated 
with the matrix in ( fTOl i is strongly connected. This implies that 
the matrix is irreducible [28|. Since the matrix is irreducible, 
weakly diagonally dominant, and has at least one row with 
strict diagonal dominance, the matrix is irreducibly diagonally 
dominant. By the Levy-Desplanques theorem, the matrix is 
non-singular ll28l . Thus, Jn is non-singular. This proves 
that the Jacobian for a connected, lossless system at a zero 
power injection solution is invertible, under the assumption of 
inductive lines. 

Although the assumptions of lossless systems and inductive 
hnes are required for the above proof, non-singularity of the 
Jacobian at a zero power injection solution generically holds 
for more general systems (e.g., lossless systems with some 
capacitive lines and lossy systems). A singular Jacobian would 
imply marginal stability at the zero power injection solution 
with multiple solutions coalescing at a bifurcation point. There 
is no reason to expect this to occur at a zero power injection 
solution. Zero power injection solutions for all IEEE power 
flow test cases f2E\ have non-singular Jacobians at zero power 
injection solutions. 

If the Jacobian of the power flow equations is non-singular 
at the zero power injection solution, the implicit function 
theorem indicates that a solution must persist for all power in- 
jections in a small ball around the zero power injection. Thus, 
there exists some voltage magnitude and angle perturbation 
AVZA5 such that 



IV. Sufficient Condition for Power Flow 
Insolvability 

A. Condition Description 

The proof in Section |lll] shows that there exists a voltage 
profile satisfying the power injection equations. We develop 
the sufficient condition for power flow insolvability by de- 
termining whether any such voltage profile could match the 
specified slack bus and PV bus voltages. No solution to the 
power flow equations exists if it is impossible to obtain a 
voltage profile that yields the specified power injections while 
also matching the specified voltage magnitudes at slack and 
PV buses. 

One way to determine if a valid voltage profile exists is 
to find the voltage profile with the lowest possible slack bus 
voltage. If the minimum possible slack bus voltage is greater 
than the specified slack bus voltage, no voltage profile will 
satisfy the power flow equations and thus the power flow 
equations are insolvable. If, conversely, the minimum slack 
bus voltage is less than the specified slack bus voltage, a power 
flow solution may exist. 

This condition thus indicates that no power flow solution 
exists when the minimum slack bus voltage obtainable while 
satisfying the power injection equations (with PV bus voltage 
magnitudes scaled proportionally) is greater than the specified 
slack bus voltage magnitude. An optimization problem with 
objective function minimizing the slack bus voltage and con- 
straints on power injections and PV bus voltage magnitudes, 
as shown in ( fTSl l. is used to evaluate this condition. 



f {V + AVZAS) = AP + jAQ (11) 

for any small AP, AQ, where V is the voltage profile for 
the zero power injection solution, AP and AQ are small 
perturbations to the active and reactive power injections, and 
/ represents the power flow equations relating the voltages 
and power injections. 

C. Scaling Up Voltages 

We complete the solution existence proof by expanding the 
small ball around the zero power injection solution to obtain a 
voltage profile that yields the originally specified power injec- 
tions. Since the power flow equations are quadratic in voltage 
magnitudes V, scaling all voltage magnitudes also scales the 
power injections. That is, scaling the voltage magnitudes in 
(fTTT l by the scalar (3 gives 

f {(3 {V + AVZA5))^ (AP + jAQ) (12) 

Choose a AP+jAQ that is in the direction of the specified 
power injections and obtain a corresponding voltage profile 
V+ AVZAS. Then increase /3 until the power injections given 
by f {P {V + AVZAS)) match the specified power injections. 
The voltage profile /3 {V + AV^ZA^) then yields the specified 
power injections. 



min Vsiack (13a) 
subject to 

n 

Pk = Vfe ^ Vi (Gjfc cos {5k - Si) + B,fc sin (4 - S^)) 

'=1 yk€{PQ,PV} (13b) 

n 

Qfe = Vfc ^ V^ {Gik sin {5k - 5,) - B^- cos (4 - 5^)) 

"ikeVQ (13c) 

14 = afcV;iack VfcePV (13d) 

where P Q is the set of PQ buses, VV is the set of PV buses, 
and Vsiack is the slack bus voltage magnitude, represents 
the specified ratio of the PV bus k and slack bus voltage 
magnitudes. The minimum achievable slack bus voltage (i.e., 
the optimal objective value of (fTsl l) is denoted as V^^^j^. 

Note that (flSl l is feasible for lossless systems and expected 
to be feasible for most practical power systems as argued 
in Section HII] The power flow insolvability condition can 
therefore be evaluated. 

The optimization problem ( fT3] l is in general non-convex 
||241 . and hence solution for a global optimum is not assured. 
A global minimum is required in order to ensure the validity of 
the sufficient condition on power flow solution non-existence. 
We therefore formulate in (Il4l the semidefinite dual of (fTjI . 
SDP algorithms can assure that we find the global solution to 
the convex dual formulation (fl4] i. 
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max ^ (AfePfc) + ^ (ifcQfc 

ke{VQ,vv} kevQ 
subject to 



(14a) 



A(A,7,/i) = 



Mslack - {^''^k + 7fcYfc) 



J2 (^feYfc + flk (Mfe - OfcMslack)) 



y (14b) 



where A/c, 7fe, and /ife are the Lagrange multipliers for active 
power (equation ( I13bl i). reactive power (equation (I13ct ). and 
PV bus voltage magnitude ratio (equation (I13db ) constraints, 
respectively, associated with bus k. The symbol ^ indicates 
that the corresponding matrix is constrained to be positive 
semidefinite. The maximum lower bound on the minimum 
achievable slack bus voltage (i.e., the optimal objective value 
of Gil) is denoted as HXcfe. 

Matrices employed in (fl4] l are defined as 



Re (n + ) Im {Y^ - Y^) 



Im(n-Yfe^) Re{Yk +Y, 



T 



Im {Yk + Y^) Re - 
Re {Y^ - n) Im + Y^) 




T 



(15) 
(16) 
(17) 



where denotes the fc*'' standard basis vector in M" and the 
matrix = e^e^Y. Notation is adopted from ll22l . 

Note that the dual formulation (fl4] l is always feasible since 
the point A,; = 0, 7,; = 0, /i^ = for all i implies A = 

y^-alack t 0. 

The semidefinite dual formulation (fl4l l provides a lower 
bound on the minimum slack bus voltage in ( fT3] l. No solution 
to the power flow equations exists if the lower bound from 
(fl4] i is greater than the specified slack bus voltage. That is, 
the condition 



V] 



lack 



>Vo 



(18) 



where Vq is the specified slack bus voltage, is a sufficient 
condition for insolvability of the power flow equations. Note 
that this formulation does not use any rank relaxations or 
enforce any requirements on matrix ranks; the solution to the 
convex problem ( fT4l i is only used as a lower bound on ( fT3] l. 

The converse condition does not necessarily hold; the power 
flow equations may not have a solution even if 



rrrmn ^ -,r 



(19) 



Thus, (fT9] l is a necessary, but not sufficient, condition for 
power flow solvability. However, satisfaction of ([19) is ex- 
pected to often predict the existence of a power flow solution. 

If the A matrix in ( I14bb has a nullspace with rank less 
than or equal to two, a solution of the power flow equa- 
tions with slack bus voltage equal to 11^™^. (and PV bus 
voltage magnitudes scaled proportionally) can be obtained 
(see II22I for further details). If a solution to the power flow 



equations with slack bus voltage equal to Vq does not exist, 
the solution with lower slack bus voltage must disappear 
as the controlled voltages increase. The disappearance of 
a solution due to increasing controlled voltages does not 
typically occur Thus, satisfaction of (fT9] l by a solution to (fl4b 
with rank (null (A)) < 2 is a strong indicator of power flow 
solution existence. 

B. Controlled Voltage Margin 

The sufficient condition (fTsT l is binary; the specified power 
flow equations either cannot have a solution or may have a 
solution. We next show how the sufficient condition can be 
interpreted to give a measure of the degree of solvability. 
We develop a measure of the distance to the power flow 
solvability boundary, which we define as the set of solvable 
power injections where all solutions may vanish under small 
perturbations. Since operating a power system far from the 
power flow solvability boundary is required to ensure stability, 
a measure of the distance to the solvability boundary is useful. 
A measure of the distance to the solvability boundary also 
indicates how close insolvable power flow equations are to 
solvability. 

We introduce a controlled voltage margin measure a for the 
distance to the power flow solvability boundary. The controlled 
voltage margin is defined as the ratio between the specified 
slack bus voltage and the lower bound on the minimum slack 
bus voltage obtained from (fl4l i. 



yrntn 
—stack 



(20) 



(7 is an upper (non-conservative) bound of the distance to 
the power flow solvability boundary. For a solvable set of 
power flow equations, we are guaranteed to be at or beyond 
the solvability boundary if the specified slack bus voltage 
decreases by the factor a. For insolvable power flow equations, 
increasing the slack bus voltage magnitude (with proportional 
increases in PV bus voltage magnitudes) by at least a factor 
of i (without changing the power injections) is required for 
solvability. 

The sufficient condition can be rewritten in terms of the 
voltage margin; ct < 1 is a sufficient condition for power flow 
insolvability. 

C. Power Injection Margin 

Another measure of the distance to the power flow solvabil- 
ity boundary is given in terms of power injections. The power 
injection margin is a measure of how large of a change in the 
power injections in a certain profile is required for the power 
injections to be on the solvability boundary. We consider the 
profile where power injections are uniformly changed at each 
bus in order to take advantage of the quadratic nature of the 
optimization problem (fTsl i in the sufficient condition. (The 
impact of non-uniform changes in power injections is future 
work discussed in Section IVIII ) The quadratic property that 
we exploit can be written as 



h{v{P + jQ))=v{v7aZky 



(21) 
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where P and Q are vectors of the active and reactive power 
injection at each bus, h is the function representing optimiza- 
tion problem (fT3] l relating the minimum slack bus voltage to 
the power injections, and 77 is a scalar 

( |2TI ) describes the linear relationship between the square 
of the voltages and the power injections. This relationship is 
evident from ( |13bt and (|13c| l: scaling all voltages by ^ scales 
the active and reactive power injections by 77. 

To develop the power injection margin, uniformly scale the 
power injections until the sufficient condition (fTSl l indicates 
that the power injections are (at least) on the solvability 
boundary. 

v{vT,::,)' ^{Vof (22) 

The power injection margin 77 corresponding to the con- 
dition in ( |22] | gives an upper, non-conservative bound of 
the distance to the solvability boundary in the direction of 
uniformly increasing power injections. For a solvable set of 
power injections, the largest proportional increase in power 
injections at each bus while potentially maintaining solvability 
is a factor of 1]. For an insolvable set of power injections, a 
proportional change of all power injections by at least r/ is 
required for a solution to be possible. 

Note that the power injection margin can be rewritten in 
terms of the voltage margin. 

V = i^f (23) 

The sufficient condition for power flow insolvability can be 
rewritten in terms of the power injection margin: < 1 is a 
sufficient condition for power flow insolvabiUty. 

The power injection margin can alternatively be calculated 
using a different optimization problem that directly maximizes 
power injections at constant power factor with fixed slack 
and PV bus voltage magnitudes. This gives the same power 
injection margins as (|22] | in all test cases. This approach lacks 
a feasibility proof and does not have any advantages over (|22T) . 

V. Numeric Example 

We next apply the sufficient condition for power flow 
solvability to the IEEE 14-bus and IEEE 118-bus systems 
Il26| using optimization codes YALMIP |29| and SeDuMi | 30]. 
The power injections are uniformly increased at each bus at 
constant power factor until the sufficient condition indicates 
that no solutions are possible. The sufficient condition results 
are compared to power flow solution attempts by a Newton- 
Raphson algorithm. 

A. IEEE 14-Bus System Results 

Results from applying the sufficient condition to the IEEE 
14-bus system are given in Table |T] The specified slack bus 
voltage is Vb = 1.0600 per unit. 

The originally specified active and reactive power injections 
are increased uniformly at each bus. The first column of Table 
H] lists the multiple by which the injections are increased. No 
power flow solutions exist after a sufficiently large increase 
(approximately 4.060 for this example). Note that the injection 



Injection Multiplier 


NR Converged 


— slack 


1.000 


Yes 


0.5261 


2.000 


Yes 


0.7440 


3.000 


Yes 


0.91 12 


4.000 


Yes 


1.0522 


4.010 


Yes 


1.0535 


4.020 


Yes 


1.0548 


4.030 


Yes 


1.0561 


4.040 


Yes 


1.0575 


4 050 


Yes 


1 0588 


4.055 


Yes 


1.0594 


4.056 


Yes 


1.0595 


4.057 


Yes 


1.0597 


4.058 


Yes 


1.0598 


4.059 


Yes 


1.0599 


4.060 


No 


1.0601 


4.061 


No 


1 .0602 


4.062 


No 


1.0603 


4.063 


No 


1.0605 


4.064 


No 


1.0606 


4.065 


No 


1.0607 


5.000 


No 


1.1764 



TABLE I 

Solvability Condition Results For IEEE 14-Bus System 



multiplier given in the first column does not change at a 
constant rate but rather focuses on the region near power flow 
solution non-existence. 

The second column indicates whether a Newton-Raphson 
solver converged to a solution at the corresponding loading. 
In order to increase the likelihood of convergence, the Newton- 
Raphson solver was initialized at each injection multiplier with 
the solution from the previous injection multiplier and a large 
number of Newton-Raphson iterations were allowed. 

The third column provides the lower bound on the minimum 
slack bus voltage in per unit obtained from ([14). In order to 
evaluate the sufficient condition for power flow insolvability at 
each injection multiplier, the value in this column is compared 
to the specified slack bus voltage of 1.06 per unit. If the value 
in the third column is greater than 1 .06, the sufficient condition 
indicates that no power flow solutions exist. These results 
show agreement between Newton-Raphson convergence and 
the sufficient condition; a power flow solution was found 
for all injection multipliers where the sufficient condition 
indicated that a solution was possible (observe that both V^i"^^. 
is just greater than 1.06 and no solution is found by the 
Newton-Raphson solver at an injection multiplier of 4.060). 

The existence of a solution for all power injections that 
satisfy (fT9] l is expected since the A matrix in ( I14bb has a 
nullspace with rank two. This need not always be the case. 
In IV-BI we investigate an example with rank (null (A)) = 4 
where no solution was found for some power injections 
even though the condition (T% indicated that a solution was 
possible. 

We next use the IEEE 14-bus system example to demon- 
strate the voltage and power injection margins. In Fig. |2] the 
voltage margin cr is plotted versus the injection multiplier 
The voltage margin decreases as power injections increase. 
The voltage margin crosses one at an injection multiplier of 
4.0595, indicating that no power flow solution can exist for 
larger power injections. Beyond this point, the voltage margin 
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IEEE 1 4-Bus Voltage Margin a vs. Injection Multiplier 



IEEE 14-Bus Continuation Trace: Bus 5 




2 3 4 

Injection Multiplier 



Fig. 2. IEEE 14-bus Voltage Mat-gin 



provides the minimum increase in the slack bus vohage (with 
corresponding proportional voltage increases at all PV buses) 
required in order for a power flow solution to possibly exist. 

In Fig. [3] we examine the power versus voltage (PV) 
curves for the high-voltage, stable solution to the IEEE 14-bus 
system. These curves, which were plotted using continuation 
techniques [311, show how the solution voltages change with 
proportional increases in power injections at all buses. The 
plots show the voltage at the arbitrarily selected PQ bus five. 
(Plotting the voltage at a PQ bus is required since voltage 
magnitudes at slack and PV buses are fixed.) The PV curve 
using the nominal slack and PV bus voltages is shown in black. 

Evaluating the optimization problem (fT4] i at an injection 
multiplier of one gives a y_ 



1.0600 _ 
0.5261 



X'fc = 0.5261. The voltage margin 
= 2.0148 per unit. Thus, no solution can exist 
if the slack bus voltage is reduced by more than a factor of 
2.0148 (with all PV bus voltages reduced proportionally). The 
grey PV curve in Fig. |3a]is obtained when the voltages are thus 
reduced. This curve shows that with these reduced voltages, 
there is the single solution is on the power flow solvability 
boundary; no solutions exist after any further increase in 
the injection multiplier Thus, the voltage margin accurately 
indicates the distance to power flow insolvability. 

The solution to the optimization problem (fT4] l also enables 
determination of the power injection margin r/. Solving (l22T i 
yields ij = (g^ff^x) = 4.0595. Thus, the power injections can 
be increased uniformly by a factor of 4.0595 until the sufficient 
condition indicates that no power flow solutions are possible. 
The black PV curve associated with the nominal voltages in 
Fig. ISalcorroborates this assertion: a power flow solution exists 
for all power injection multipliers less than 4.0595, but no 
solution exists beyond this power injection multiplier 

The voltage and power injection margins can also be used to 
investigate insolvable power injections. Assume that we would 
like to consider operation at a power injection multiplier equal 
to five. Evaluating the optimization problem (fT4l i at a power 
injection multiplier of five gives Vj^ia^k = 1-1764. Note that 
{22i implies that knowledge of Hs^acfe ^ power injection 
multiplier of one allows the direct calculation V^^^. at a 
power injection multiplier of five: 




2 3 4 

Injection Multiplier 



1 

(a) Nominal and Low Slack Bus Voltage PV Curves 
IEEE 14-Bus Continuation Trace: Bus 5 




2 3 4 5 

Injection Multiplier 

(b) Nominal and High Slack Bus Voltage PV Cuires 
Fig. 3. IEEE 14-Bus System PV Curves 



V" 



^rmn 
slack I 



I Inj Mult — 5 I Inj Mult- 

= 0.5261 ■ %/5 = 1.1764 perunit 



(24) 



The voltage margin at a power injection multiplier of five is 
c — 1 ^'^L = 0.9011. (T < 1 indicates that there is no solution 

1.1764 

at a power injection multiplier of five. To potentially achieve 
a power flow solution, the slack bus voltage must increase 
by at least a factor of Q-gfjjj- = 1.1098 (with corresponding 
proportional increases in all PV bus voltages). The grey PV 
curve in Fig. |3b]has the voltages thus increased. Observe that 
increasing the voltages in this way allows for a power flow 
solution on the power flow solvability boundary for a power 
injection multiplier of five. 

The power injection margin rj can also be calculated at a 
power injection multiplier of five using 



Vo 



-!-s!acfe llnj Mult=5 



1.0600 y 

1.1764 / 



: 0.8119 



(25) 



77 < 1 implies that no solution exists at a power injection 
multiplier of five. The power injection margin also indicates 
that no solution can exist for power injection multipliers 
greater than 0.8119 • 5 = 4.0595. This corresponds to the 
"nose" point of the black (nominal) PV curve in Fig. |3b] 
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Injection Multiplier 


NR Converged 


ymin 
— slack 


1.00 


Yes 


0.5724 


1.50 


Yes 


0.7010 


2.00 


Yes 


0.8095 


2 50 


Yes 


9050 


3.00 


Yes 


0.9914 


3.15 


Yes 


1.0159 


3.16 


Yes 


1.0175 


3.17 


Yes 


1.0191 


3.18 


Yes 


1.0207 


3.19 


No 


1 .0223 


3.20 


No 


1.0239 


3 21 


No 


1 0255 


3.22 


No 


1.0271 


3.23 


No 


1.0287 


3.24 


No 


1.0303 


3.25 


No 


1.0319 


3.26 


No 


1.0335 


3.27 


No 


1.0351 


3.28 


No 


1.0366 


3.29 


No 


1.0382 


4.00 


No 


1.1448 



TABLE II 

INSOLVABILITY CONDITION RESULTS FOR IEEE 1 1 8-BUS SYSTEM 



B. IEEE 118-Bus System Results 

Results from applying the sufficient condition to the IEEE 
1 18-bus system are given in Table HIl The data are arranged in 
the same manner as in Table J] The specified slack bus voltage 
is Vb = 1.0350 per unit. 

The results can be categorized intro three regions: small 
power injections where the sufficient condition indicates that 
a solution is possible and a solution is indeed found using 
a Newton-Raphson solver, larger power injections where the 
sufficient condition indicates that a solution is possible but no 
solutions are found, and yet larger power injections where the 
sufficient condition indicates that no solutions are possible and 
no solutions are found. 

The rank of the nullspace of the A matrix in (I14bl i for the 
IEEE 118-bus system was four Therefore, the expectation that 
satisfaction of ( fT9l l will result in power flow solvability may 
not hold. Correspondingly, these results emphasize the fact 
that ( fTSl l is a sufficient condition for power flow insolvability. 
Specifically, a power flow solution was not found for injection 
multipliers greater than 3.18, even though V^siack ^^^^ '^^an 
the specified slack bus voltage until an injection multiplier of 
3.27. A continuation power flow indicates that the high-voltage 
solution bifurcates at a power injection multiplier of 3.185, so 
it is likely that no solutions exist after this point. No solutions 
are found at injection multipliers larger than 3.27 where the 
sufficient condition indicates that no solutions are possible. 

We next use the IEEE 118-bus system example to demon- 
strate the voltage and power injection margins. In Fig. |4] the 
voltage margin cr is plotted versus the injection multiplier 
Similar to Fig. |2] the voltage margin decreases as power 
injections increase. The voltage margin crosses one at an 
injection multiplier of 3.2695, indicating that no power flow 
solution can exist for larger power injections. For larger power 
injections, the voltage margin provides the minimum increase 
in the slack bus voltage (with corresponding proportional 



IEEE 1 1 8-Bus Voltage Margin g vs. Injection Multiplier 




O' ' ' ' 

12 3 4 

Injection Multiplier 
Fig. 4. IEEE 118-bus Voltage Margin 

IEEE 118-Bus Continuation Trace: Bus 44 




0.5 1 1.5 2 2.5 3 3.5 
Injection IVIultiplier 

(a) Nominal and Low Slack Bus Voltage PV Curves 



IEEE 118-Bus Continuation Trace: Bus 44 

1.2^ - - 




12 3 4 

Injection Multiplier 



(b) Nominal and High Slack Bus Voltage PV Curves 
Fig. 5. IEEE 118-Bus System PV Curves 



voltage increases at all PV buses) that is required in order 
for a power flow solution to possibly exist. 

In Fig. |5] we examine the PV curves for the high-voltage, 
stable solution to the IEEE 118-bus system. The plots show 
the voltage at the arbitrarily selected PQ bus 44 (plotting the 
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voltage at a PQ bus is required since the voltages at slack and 
PV buses are fixed). The PV curve using the nominal slack 
and PV bus voltages is shown in black. 

Evaluating the optimization problem (fl4] l associated with 
the sufficient condition at an injection multiplier of one gives 
0.5724. The voltage margin is a 



-irnmi 
<^ jLslack 



0.5724 



1.8082. Thus, no solution can exist if the slack bus voltage 
is reduced by more than a factor of 1.8082 (with all PV bus 
voltages reduced proportionally). The grey PV curve in Fig. |5a] 
is obtained when the voltages are thus reduced. Although no 
solutions exist for the grey PV curve at injection multipliers 
larger than one, there are also injection multipliers slightly less 
than one for which no solutions are found with continuation 
techniques. This reinforces the fact that the voltage margin is 
an upper bound on the distance to the solvability boundary. 

The solution to the optimization problem (fl4] i also enables 
determination of the power injection margin 77. Solving (l22T l 
yields 77 = (g^§f§f )" = 3.2695. Thus, the power injections can 
be increased uniformly by a factor of 3.2695 until the sufficient 
condition indicates that no power flow solutions are possible. 
This is also an upper bound on the distance to the solvability 
boundary: as indicated by the sufficient condition, the black 
PV curve associated with the nominal voltages in Fig. ISalhas 
no solutions for power injection multipliers larger than 3.2695, 
but also appears to have no solutions for some values of power 
injection multipliers below 3.2695. (It is possible, but unlikely, 
that a PV curve associated with a different solution may exist 
at injection multipliers between the "nose" of the PV curve 
associated with the high-voltage solution at 3.1840 and the 
value of 3.2695 from the sufficient condition.) 

The voltage and power injection margins can also be used to 
investigate insolvable power injections. Assume that we would 
like to consider operation at a power injection multiplier equal 
to four Evaluating the optimization problem (fT4l i at a power 
injection multiplier of four gives VlXcfc 1-1448. Note that 
( |22] | implies that knowledge of V^i"^.^. at a power injection 
multiplier of one allows the direct calculation V^^j^ at a 
power injection multiplier of four: 



xrrmn 

-!^s;acfc I inj Mult=4 



r- Y-min I 
V'/ -!^s;acfc|lnjMult = l 

0.5724 •\/4 = 1.1448 per unit (26) 



The voltage margin at a power injection multiplier of four is 
cr = Yaffil = 0.9041. (T < 1 indicates that there is no solution 
at a power injection multiplier of four To potentially achieve 
a power flow solution, the slack bus voltage must increase 
by at least a factor of g^j^jgg = 1.1061 (with corresponding 
proportional increases in all PV bus voltages). The grey PV 
curve in Fig. |5b] has the voltages thus increased. Since no 
solutions are evident from the PV curve at an injection multi- 
plier of four, it appears that this is not a large enough voltage 
increase to obtain solvability. This is a result of the fact that 
we use a sufficient condition for power flow insolvability to 
calculate the voltage margin; failing to satisfy the the sufficient 
condition for power flow insolvability does not ensure the 
existence of a solution. 



The power injection margin rj can also be calculated at a 
power injection multiplier of four using 



, —siacfe I i,ijMult=4/ 



1.0350 



1.1448 / 



0.8174 



(27) 



77 < 1 implies that no solution exists at a power injection 
multiplier of four The power injection multiplier also indicates 
that no solution can exist for power injection multipliers 
greater than 0.8174 • 4 = 3.2695. The "nose" point of the 
black (nominal) PV curve in Fig. |3b] is slightly lower than 
this upper limit on power flow solvability. 

VI. Solutions with Non-Zero Duality Gap 

If the A matrix in ( |14b| l has nullspace with rank greater than 
two, a non-zero duality gap exists between the primal problem 
( fT3] l and the dual problem (fl4] l ll22l . The dual problem will 



give a F™acfe that is less than the V™™^. of the primal problem. 
If H'Xcfe <Vo< VXcfe' condition ^ may hold without the 



existence of a power flow solution. Since directly calculating 
^s'lack difficult due to the non-convexity of the power flow 
equations in ( fT3]) . it is hard to determine whether a solution 
exists for cases where (fT9t holds and rank (null (A)) > 2. The 
IEEE 118-bus system in Section [V-BI provides an example of 
such a case. 

As the controlled voltages are decreased, the solutions to 
the power flow equations come together in bifurcations. We 
conjecture that the rank of the nullspace of the A matrix 
depends on where these solutions bifurcate. The minimum 
slack bus voltage of the primal optimization (fljT l is obtained 
when the last solutions bifurcate. If the minimum slack bus 
voltage is obtained when two solutions bifurcate far from 
the point of bifurcation of any previously existing solutions, 
we conjecture that rank (null (A)) < 2. If more than two 
solutions bifurcate at or near the same minimum slack bus 
voltage (e.g., the final two pairs of solutions bifurcate with 
the same or similar values of Vy^")^), we conjecture that 
rank (null (A)) > 2. 

The system diagram in Fig.|6]provides an example with two 
pairs of solutions bifurcating at the same value of minimum 
slack bus voltage, resulting in rank (rank (A)) — 4. Bus one 
is a PV bus with voltage magnitude 1.0 per unit and generation 
of 0.5 per unit, bus two is a slack bus with voltage 1.0Z0°, 
and bus three is a PQ bus with load 1.0 + jO.25 per unit. 

Fig. I2] shows the slack bus reactive power injection Q2 for 
all four existing solutions as the controlled voltages at buses 



|vj = 1.0 



V =1.0Z.o" 




= 1.0+j0.25 



R=0.5 



Fig. 6. Three-Bus System 
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Slack Bus Reactive Power Injection vs. Controlled Voltages 




0.4 0.5 0.6 0.7 0.8 0.9 1 
Fig. 7. Solutions for Three-Bus System 

one and two are uniformly decreased. When the controlled 
voltage magnitudes reach 0.3582 per unit, both pairs of 
solutions bifurcate. Note that both pairs of solutions bifurcate 
at different values of slack bus reactive power injection. That 
is, even though both solution pairs disappear at the same value 
of controlled voltages, all four solutions do not come together 
to the same point. 

For this example, V^"/™^. from optimization (fT4l i finds the 
correct minimum slack bus voltage of 0.3582 per unit. The 
bifurcation points for both solution pairs in Fig. Q have this 
voltage. Thus, the optimization problem ( fT4] i has no bias 
toward choosing either of these two solutions. Indeed, since all 
solutions must be connected in the convex dual formulation, 
any matrices in the dual feasible space that lie between these 
solutions, including matrices that have nullspaces with high 
ranks, are possible solutions to (fl4l l. Thus, the optimization 
problem (fl4] i is not expected to obtain rank (null (A)) < 2. 
We conjecture that similar explanations apply to other systems 
with rank (null (A)) > 2, such as the IEEE 118-bus system 
used in Section [V-BI 

Using optimization codes YALMIP fllgl and SeDuMi f30l . 
we obtained a solution to (fT4l i for the three-bus system 
with rank (null (A)) = 4. For this system, adding a small 
bias in the objective function to favor larger or smaller 
slack bus reactive power injection resulted in solutions with 
rank (null (A)) = 2. Appropriate biases to obtain solutions 
with rank (null (A)) = 2 have not yet been found for other 
systems (e.g., the IEEE 118-bus system used in Section [V-BI ). 

The sufficient condition (fTsT l does not generally give a 
solution to the power flow equations with slack bus voltage 
equal to the specified value of Vb; traditional techniques are 
required to actually find a potentially existing solution. In an 
attempt to actually find power flow solutions, a modification 
to the objective function of (fT3T l was attempted. Specifically, 
the objective function (I13ab was replaced with (|28] |. 

min (Vl,,- 1/0^)2 (28) 
This modified objective function minimizes the squared 



difference in squared slack bus voltage magnitude (squared 
voltage magnitudes are necessary to easily form the dual 
problem). The power flow equations have no solution when 
the optimal objective value is greater than zero. 

Since the constraints ( I13bb . ( I13cl i. (I13dl i are unchanged, the 
existence proof given in Section |lll] applies, so this modified 
condition can be evaluated. 

In the same manner as the original sufficient condition, the 
dual of the optimization problem with the modified objective 
function ( l28l l provides a lower bound on the primal objective 
function. If the dual objective value is greater than zero, no 
solutions exist for the power flow equations. If, additionally, 
the A matrix has a nullspace with rank less than or equal to 
two, the solution to the power flow equations can be obtained 

Unfortunately, the optimization problem with this modified 
objective function generally failed to provide A matrices that 
had nullspaces with appropriate rank on all IEEE power flow 
test systems 1261 . Also, the modified formulation often had 
numeric problems. This modified objective function (|28) was 
thus not practically useful. Future work includes investigating 
other modifications that allow for determining the solution to 
the power flow equations in addition to providing a sufficient 
condition of solution non-existence. 

VII. Conclusion and Future Work 

We have presented a sufficient condition for identifying 
insolvability of the power flow equations. This sufficient 
condition requires the evaluation of an optimization problem. 
We have proven that this optimization problem is feasible 
for lossless power systems and argued that practical power 
systems should also yield a feasible optimization problem. 
In order to quantify the degree of solvability, we developed 
controlled voltage and power injection margins from the 
sufficient condition that provide upper bounds on the distance 
to the power flow solvability boundary. Finally, we applied 
the sufficient conditions, voltage margin, and power injection 
margin to the IEEE 14-bus and 118-bus example systems. 
The IEEE 118-bus system exemplified the fact that this is a 
sufficient condition for power flow insolvability; it is possible 
to find cases that appear to have no solution even though the 
sufficient condition for insolvability is not satisfied. However, 
the majority of power systems we investigated yielded results 
similar to the IEEE 14-bus system where a power flow solution 
was found with a Newton-Raphson algorithm up to the point 
identified by the sufficient condition as insolvable. 

There are several open questions associated with this work 
that we intend to pursue. We will investigate how the voltage 
and power injection margins change with non-uniform power 
injection scaling. For instance, what happens if the power 
injections are scaled at a subset of the buses, or if the 
power injections are scaled with non-constant power factor? 
The sensitivity of V™™^ in ( fT3] l to the active and reactive 
power flow constraints (I13bt and (I13cb (as determined by the 
Lagrange multipliers A and 7, respectively) provides insight 
for small changes, but numerical experience indicates that 
these sensitivities are only valid for very small perturbations. 



12 



University of Wisconsin-Madison Department of Electrical and Computer Engineering. Technical Report ECE-12-01. 



The impact of large power injection changes requires further 
investigation. 

Another area open to further study regards computational 
issues. The size of positive semidefinite constraint in the 
dual optimization problem ( I14bl i scales as (2n)^, where n 
is the number of buses in the system. This scaling prohibits 
directly solving (fT4l l for large systems. Future work includes 
researching whether the structure and sparsity of the power 
flow equations can be exploited to solve this optimization 
problem quickly and efficiently for large systems. 

Finally, we intend to further investigate our conjectured ex- 
planation of solutions with non-zero duality gap. Do multiple 
pairs of solutions bifurcate near the minimum slack bus voltage 
for other systems? What system topologies yield the possibility 
of having rank (null (A)) > 2? 
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